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Abstract 

o 

D ' The global response to weak time periodic forces of an array of noisy, coupled nonlinear systems 



might show a nonmonotonic dependence on the number of units in the array. This effect has been 
termed system size stochastic resonance. In this paper, we analyze the nonmonotonic dependence 
on the system size of the signal-to-noise ratio of a collective variable characterizing a finite array 
of one-dimensional globally coupled bistable elements. By contrast with the conventional non- 



^ . monotonic dependence with the strength of the noise (stochastic resonance), system size stochastic 

o ; 

CJ , resonance is found to be restricted to some regions in parameter space. 
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I. INTRODUCTION 



In some nonlinear systems, noise can be used to improve the response of the system to 
a weak external signal. In the last two decades, this rather counterintuitive effect known 
as Stochastic Resonance (SR), has been massively studied, both experimentally and the- 
orefcany. due to ,ts potential interest to a variety of areas of scenee and technology 
[2|. A large majority of the experimental as well as theoretical studies have dealt mainly 
with the response of single systems to weak external forcings. Nonetheless, the response 



of complex systems 
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ie arrays composed of many coupled units has also been considered 
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12|. 



The study of complex systems opens new possibilities with respect to those observed in 
single unit systems. Besides the usual noise strength, the variation of other system parame- 
ters can give rise to enhanced, nonmonotonic behaviors of some magnitudes characterizing 
the system response to a weak externally applied time periodic force. Indeed, nonmonotonic 
behaviors of pertinent quantifiers with the coupl ing strength parameter of the array have 

nnnnflnnnn 

been reported in the literature [3|, 14, l5|, la, 17|, l8|, l9|, [ly. 111! ]. 

The noise enhanced nonmonotonic response of driven systems with the number of units 
has been termed system size stochastic resonance (SSSR). It has been studied in several con- 
texts. Several authors have considered the dependence of system properties on the number 

nnn 

of elementary units in complex arrays of bistable elements [1^, Il3|, [IJ]. The dependence 
on the number of ion channels of the ion concentrations along cell membranes in biological 



17|, the authors analyze the 



Hodgkin-Huxley type models has been studied in In 
optimization of neuron firing events in coupled excitable neuron units as the number of 
coupled elements is varied. The phenomenon of SSSR has also been discussed within the 
context of modelling opinion formation in social collectivities fisl . 

In this paper, we address the phenomenon of SSSR for arrays of globally coupled noisy 
bistable units. We consider that the collective response of the whole complex system is 
well represented by a single random variable. Its dynamics is clearly constructed from the 
dynamics of the constituent elements of the array but it shows collective emergent properties. 
In Section [IT] we introduce a model in terms of the individual units and the interactions among 
them, and we define the collective variable and the quantifier that we will use to characterize 
the phenomenon. Next, in Section UTTl we discuss the possibility of a reduced description of 
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the collective variable dynamics in terms of an effective Langevin equation, as put forward 
by Pikovsky et al. in 12]. In Section HVl we address the problem with numerical simulation 
tools and discuss the results obtained. The paper ends with some conclusions. 



II. THE MODEL 

The model considered in this work is the same as the one introduced by Desai and 
Zwanzig in their influential paper [l9|, except that we add a time periodic force to the 
1/angevin dynamics of each degree of freedom. This model is also the one considered in 
1^ . The model consists of a set of N identical bistable units, each of them characterized 
by a variable Xi{t) {i = 1, . . . ,N) satisfying stochastic evolution equations (in dimensionless 
form) of the type 

6 ^ 

Xi = Xi-x^,+- Y.{x, - Xi) + V2DUt) + F{t), (1) 
i=i 

where ^ is a coupling parameter and the term ^j(t) represents a white noise with zero average 
and {C,i(t)C,j{s)) = 6ij6{t — s). The external driving force is periodic in t, F(t) = F(t + T). 
An alternative formulation of the dynamics is in terms of the Fokker-Planck equation for 
the joint probability density f^ixi, . . . ,XN,t), 

where ?7 is a time dependent potential energy relief. 



AT n N N 



Q 

U{Xi,...,XN,t) = Y^^i^i) + J^J2Y.(^i - ^i, 

i=l i=l j=l 



2 



N 

+F{t)Y,x, (3) 



with 
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We are interested in the properties of a collective variable, S{t), defined as 

1 ^ 

S{t) = -Y.^,{t), (5) 
i=i 

characterizing the chain as a whole. Even though the set is a non-stationary A^- 

dimensional Markovian process, S{t) is not, in general, a 1-dimensional Markovian process. 



Several magnitudes showing enhanced nonmonotonic behaviors as a system parameter is 
varied have been used as SR quantifiers Ij]. In this work, we will use the signal-to- noise 
ratio of the collective variable i^out as the relevant quantifier. In order to evaluate it, we 
need the knowledge of the one-time correlation function of the collective variable defined as 

C{r) = U^ dt{S{t)S{t + T))^, (6) 
1 Jo 

The notation (. . .) indicates an average over the noise realizations and the subindex oo 
indicates the long time limit of the noise average, i. e., its value after waiting for t long 

10|, it is convenient 



enough that transients have died out. As shown in our previous work 
to split C(r) into two parts, 

C(r) =Ccoh(r) + ancoh(r). (7) 
The coherent part, Ccoh(T), given by 

Ccoh(r) = - / dt {S{t))^{S{t + r))oo, (8) 



TJo 

is periodic in r with the period of the driving force. The incoherent part, Cincoh(T) arising 
from the fiuctuations of the output S{t) around its average value, decays to zero as r 
increases. 

The quantifier i?out is defined as 

^^^^ Cincoh(^) Cmcoh(^) 

where fl is the fundamental frequency of the driving force F{t), C'coh(^) is the corresponding 
Fourier coefficient in the Fourier series expansion of Ccoh(T), and Cincoh(^) is the Fourier 
transform at frequency Q of Cincoh(T)- 



III. APPROXIMATE DYNAMICS FOR S{t) 

The dynamics of the collective variable follows from its definition and the dynamics of 
the individual degrees of freedom. But, as S{t) is not Markovian in general, there is no 
reason why it should satisfy a 1-dimensional closed Langevin equation with effective drift 
and diffusion coefficients. For the same reason, there is no guarantee that its associated 
probability density P{s,t) would satisfy a Fokker-Planck equation. On the other hand, 
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having an approximate description in terms of a 1- dimensional stochastic process, would 
certainly be a helpful tool for the understanding of the dynamical behaviors. It is then not 
surprising that several approximate descriptions in a reduced 1-dimensional space for the 



collective variable have been proposed in the literature |12l . |20[ |. In 12|, Pikovsky et al. 
used the Gaussian truncation of an infinite hierarchy of equations for the cumulants and a 
slaving principle to construct an effective 1-dimensional Langevin equation. In the absence 
of external driving, its explicit form is 



2D 



S = aS-bS'' + ^—xit) 
with {x{'t)x{^)) = ^{t ~ ■s) and the coefficients a and b given by 



(10) 



a 
b 



l + 0.5( 



1) - 0.5 
4a 



12D: 



(11) 



2-e + ^J{2 + ey - 2ad 

Unfortunately, the coefficient b does not exist for all ranges of the 9 and D values. In Fig. 
([1]) the continuous line represents (2 + 0)^ — 24D = containing the function appearing in 
the definition of b. This line separates regions where the effective Langevin equation in Eq. 
( ITOl) either exists or not. Below that line, the b coefficient is real and the Langevin equation, 
Eq. ( ITOl) . exists, while for points above the line b is complex, and therefore, the effective 
Langevin equation is not valid. 

In general, i?out can not be evaluated analytically. But if we accept that S{t) satisfies 
a Langevin equation like Eq. (fTOj) with an extra term corresponding to the driving force 
added to its right hand side and this extra term is very weak, then a linear response theory 
might provide a valid approximation to analyze the system behavior. The problem is then 



similar to that of SR in a single bista 



obtain analytical approximations 
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^epotential and we can apply well known arguments to 



2l| for _Rout within the 



imits of linear response theory. 



12|. They find SSSR for weak 



This is the strategy followed by Pikovsky and coworkers in 
amplitude driving forces and some ranges of noise strength and coupling parameter values. 
Due to the difficulties with Eq. fllOl) discussed above, it is not clear whether the phenomenon 
exists for all parameter values or it is restricted to some particular regions of the parameter 
space. It is also interesting to relate the possible SSSR to the stochastic resonant effects in 
coupled arrays reported in other works 



10 



111 ] when the noise or the coupling constant 



are varied. 
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FIG. 1: The solid line represents (2 + 9)'^ — 241? = 0. For points above it, the coefficient 6 is a 
complex number and the effective Langevin equation ()10p is not valid. The dotted line joining 
black dots marks a numerically obtained transition line such that, above this line, the equilibrium 
distribution of the global variable is always monomodal and below it P°^{s) is multimodal (see the 
final paragraph of Section [TVl) . The data depicted were obtained for a system with = 30 units, 
but those values are very independent of the system size as long as it remains finite. 

IV. NUMERICAL RESULTS 

The lack of a reliable effective Langevin approximation valid for all regions of parameter 
space forces us to consider numerical simulations. Following the procedure indicated in 



22l |. we have numerically solved the Langevin equations for Xi{t) in Eq. ([T]) for very many 
noise realizations. The numerical solution of the A^-dimensional process is used to collect 
information about S{t) and construct the magnitudes that we need. In particular we will 
estimate numerically the long-time limit of the first two cumulant moments, the correlation 
function with its coherent and incoherent parts and histograms to estimate the probability 
density. The estimation of i?out is then a matter of numerically performing the quadratures 
indicated in its definition (see Eq. Qj). 

Let us consider a rectangular driving force: F{t) = A {F{t) = —A) if t G [nT/2, {n + 
l)T/2) with n even (odd). In what follows, we will always take the amplitude A = 0.05 and 
the fundamental frequency Q = 0.05. When such a force is applied to a single isolated noisy 
bistable, the response of the system is very well described by a linear response function, as 
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the dynamical effects arising from the distortion of the bistable potential induced by the 
driving are within the limits of a small perturbation theory. 



A. Strong coupling 



First we study the case of parameter values within the range considered in [12]. In Fig. 
([2]), the black circles correspond to the numerically obtained values of i?out as the number 
of particles N of the array is varied and 6 = 5.5, D = 1. The existence of SSSR is evident, 



in consonance with the results reported in 



12| . It should also be pointed out that the 



values of i^out are rather small, even for arrays with a size where i?out reaches its peak value. 
To understand why, we have analyzed the time behavior of the first two cumulants of the 
collective variable as well as the incoherent part of the fluctuations about its average. In 
Fig. (jnb) we depict the time evolution of the first two cumulant moments for the parameter 
values 9 = 5.5, D = 1. Even though the amplitude of the average output is about five times 
bigger than that of the driving force, the output signal is very noisy with a second cumulant 
much larger than the noise strength {2D)/30 ~ 0.066 associated to the effective Langevin 
equation, Eq. ffTOl) . The incoherent fluctuations of the collective variable are depicted in Fig. 
(jib). They are large and long-lasting. Consequently, the numerator in the fraction defining 
i?out (see Eq. (Q) should not be very large, while its denominator has a substantial value. 
The smallness of the Rout values is not surprising. 



The results of Pikovsky et al. [12(] rely on the validity of the linear response approximation. 
Within this limit, one assumes that the incoherent part of the output correlation function 
in a driven system can be safely approximated by the equilibrium correlation function of 
the same system in the absence of driving force. It is clear from Fig. ({Ud) that such an 
approximation is not too bad. Then, for the parameter values 6 = 5.5, D = 1, the linear 
response function can be used and as this function has a nonmonotonic behavior with A^, 
one finds SSSR. 

Let us now analyze what happens if we increase the noise strength to D = 2.5, while 
keeping all the other parameter values fixed. The black squares in Fig. ([2]) indicate that 
Ront increases monotonically with the size of the system. Thus, there is no SSSR for this 
noise strength even though the coupling parameter is still strong. Note that for these noise 



strength and coupling parameter values, the effective Langevin equation in Eq. (ITOI) does not 
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FIG. 2: The signal-to-noise ratio of the cohective output i?out as a function of iV for two different 
values of the noise intensity D. In both cases 6 = 5.5. The driving amplitude is A = 0.05 and the 
fundamental frequency = 0.05. 

exist. On the other hand, the numerical results depicted in Figs. ([3^) and (jUi) show that the 
second cumulant is reduced with respect to its value for D = 1, and that the incoherent part 
of the fluctuations in the driven system matches well the equilibrium fluctuations. Taking 
also into account the smallness of the driving amplitude, a linear response approximation 
should still provide a reliable description of the response of the system for large noises and 
strong couplings, even though SSSR does not exist in that region. 

It is worth noting that a nonmonotonic behavior of -Rout with the noise strength D, typical 
of stochastic resonance, does exist for strong coupling. An example is depicted in Fig. ([5]) for 
a system of = 30 elements with coupling strength 6 = 5.5, driven by a weak rectangular 
force with A = 0.05 and Q = 0.05. As expected, the values of i?out are quite small. 

B. Weak coupling 

Let us now consider cases with a much weaker coupling parameter. We will take ^ = 0.5. 
We will consider two values of the noise strength: D = 1 and D = 0.2. In Fig. ([6]) we depict 
the behavior of i^out with the system size. Clearly, there is no SSSR as the value of i?out 
increases with A^ for both noise values. 

The behavior of the first two cumulant moments of the collective variable for a driven 
system with A^ = 30 units and coupling parameter value 9 = 0.5 is depicted in Fig. ([7]) 



• • D=1 

■ ■ D=2.5 
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FIG. 3: (Color online) Time evolution of the first two cumulant moments of S{t) in a system with 
= 30 strongly coupled {0 = 5.5) units driven by a rectangular force with A = 0.05 and 0, = 0.05, 
for D = 2.5 (a) and L> = 1 (b). 

for two values of the noise strength: D = 0.2 in the lower panel, and D = 1 in the upper 
one. It is to be noted the large increase of the output amplitude relative to the weak 
driving amplitude for the D = 0.2 case, while in the higher noise situation, the amplitude 
amplification is small. The second cumulant plots indicate that the output noise level is also 
very much increased in the low D case, while it is kept very small in for the higher noise. 

The time dependence of the incoherent part of the correlation function of the collective 
variable for the driven system and the equilibrium correlation function of the same variable 
are shown in Fig. ([8]), for = 30, 6 = 0.5 A = 0.05, Q = 0.05 and two noise values. 
For the higher noise value, D = 1, the incoherent fluctuations in the driven system are 
practically identical to the equilibrium correlation function. This feature indicates that a 
linear response theory should then be an excellent approximation to analyze the response of 
the system to weak driving forces, when the noise strength is large even in a weak coupling 
situation. On the other hand, for the low noise strength case D = 0.2, the situation is 
different. The time behavior of Cincoh(^) is very different from the time behavior of Ceq(t). 
The correlations decay much faster in the driven system than in the absence of driving and 
from a smaller initial value. This indicates that the driving field, albeit weak, introduces 
a very strong distortion of the fluctuation dynamics of the collective variable during each 
cycle of the external driving, thus rendering inadequate the ideas behind a linear response 
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200 



FIG. 4: (Color online) Time evolution of the incoherent part of the correlation function (black 
line) and the equilibrium time correlation function (red line) of the collective variable in a system 
with = 30 strongly coupled { = 5.5) units for D = 2.5 (a), and D = 1 (b). The driving force 
is rectangular with A = 0.05 and = 0.05. 
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FIG. 5: Non-monotonic behavior of i?out with the noise strength D for a system with = 30 
units with global coupling strength 9 = 5.5. The driving force is rectangular with A = 0.05 and 
n = 0.05. 



approximation. This fact has been indeed analyzed previously by us in connection with SR 
in arrays of interacting particles 



23 



2J] . We have referred to this regime as a non-linear SR 



regime, characterized by the control of the fluctuations by the external driving. 
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FIG. 6: Monotonic behavior of i?out with the number of units A^. The data corresponds to ^ = 0.5 
and two different values of D. The driving force is rectangular with A = 0.05 and 0, = 0.05. 
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FIG. 7: (Color online) Time evolution of the first two cumulant moments of S{t) in a system with 
= 30 weakly coupled (0 = 0.5) units driven by a rectangular force with A = 0.05 and 0, = 0.05, 
for D = 1 (a) and D = 0.2 (b) 

In Fig. iQ we depict the nonmonotonic behavior of -Rout with D for an array of weakly 
coupled {6 = 0.5) = 30 units, driven by a rectangular force with A = 0.05 and Q = 0.05. 
The values of i?out are substantially larger than the ones obtained for the stronger coupling 
constant (compare with Fig. [5]). 

To help rationalize the results above, it is useful to consider the equilibrium distribution 
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FIG. 8: (Color online) Time evolution of the incoherent part of the correlation function (black 
line) and the equilibrium time correlation function (red line) of the collective variable in a system 
with = 30 weakly coupled {9 = 0.5) units for D = 1 (a), and D = 0.2 (b). The driving force is 
rectangular with A = 0.05 and = 0.05. 

function of the collective variable P^^{s) = {6{s — S(t))). This function can be estimated 
by constructing histograms from the long time results of the numerical simulations of Eq. 
([1]) with F{t) = 0. The data reveals the existence of a transition line separating different 
regions in the D-6 plane. The line joining the black circles in Fig. ([1]) corresponds to the 
transition line for a system with = 30 particles. The location of the line for systems 
with other number of particles do not differ substantially from the one plotted here. For 
parameter values corresponding to points above the line, the equilibrium distribution is 
always monomodal and centered around s = 0, while for points below the line, it is bimodal 
with one minimum at s = and two maxima at values ±so. The location of the minima 
depends on the values of the system parameters. Thus for a fixed N and 6, as the noise value 
is increased, the equilibrium distribution of the system would change shape from bimodal 
to monomodal. The nonmonotonic behavior of -Rout with D (SR) observed in both Figs. 
(for strong coupling) and (Q (for weak coupling) is associated the change in the shape of 
the equilibrium distribution function. By contrast, SSSR is restricted to a reduced region 
of parameter space involving rather high values of the coupling strength and relatively low 
values of the noise strength without involving a change in shape of P^^{s). SSSR is then 
connected to conditions such that a bimodal equilibrium distribution is only slightly distorted 
by the driving field. Note that even though a linear response approach might yield a reliable 
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FIG. 9: Non-monotonic behavior of i?out with D for = 30 weakly coupled (9 = 0.5) units. The 
driving force is rectangular with A = 0.05 and 0, = 0.05. 

description of the dynamics, as it happens for points above the transition line, the existence 
of SSSR is not guaranteed. In this sense, SSSR does not seem to be such a general effect as 
the usual SR nonmonotonic behavior with the noise strength. 

V. CONCLUSIONS 

In this work, we have analyzed the response of globally coupled arrays of noisy bistable 
systems driven by weak time periodic forces. We focus on the dynamics of a 1-dimensional 
collective variable characterizing the array as a whole. The lack of an adequate Langevin 
equation describing the reduced dynamics of the collective variable, prompts us to use nu- 
merical simulations involving all the individual dynamics. 

Our results indicate that SSSR exists only in limited regions of parameter space: for 
relatively low noise values and strongly coupled arrays. By looking at the behavior of the 
fluctuations in those regions, we find that the system response is adequately well described by 
a linear response approximation. Furthermore, the reduced dynamics can be well described 
by an effective Langevin equation. 

Even for strongly coupled arrays, SSSR does not exist when the noise values are large. 
The effective Langevin equation proposed in [l2| ceases to be valid. On the other hand, 
the equilibrium probability of the collective variable is, in this region of parameter values, 
monomodal. Thus, the effect of a weak applied external force on the collective dynamics 
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should still well described by a linear response function which has a monotonic behavior 
with N. 

For weak coupling, the SSSR disappears. Regardless of the noise strength, i?out mono- 
tonically increases with the number of particles in the array. For large noise values, P'^^{s) 
is monomodal. Thus, a linear response approximation should be a reliable description of the 
system response. On the other hand, for low noise values, the influence of the applied force 
on the fluctuations dynamics is so intense that the system operates in a non-linear regime. 

Our results demonstrate that in the strong coupling regime, there is a typical SR behavior 
with the noise strength D, when the system size N and the coupling constant 6 are kept 
fixed. The i?out values are, nonetheless, very small. This is in sharp contrast with the large 
enhancement of the SR effects observed in the weak coupling regime. 

The construction of effective reduced 1-dimensional dynamics for the relevant variable 

" n 



is very useful. The approaches followed in 



120( 1 are unfortunately limited to restricted 



regions of parameter space. We are currently working on obtaining approximate Langevin 
equations for the collective variable in driven systems valid for all regions in parameter space. 



along the lines initiated in 251]. 
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